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Herein, we demonstrate the control of magnetotactic bacteria through the appUcation of magnetic field 
gradients with real-time visualization. We accomplish this control by integrating a pair of macroscale 
Helmholtz coils and lithographically fabricated nanoscale islands composed of permalloy (Ni8oFe2o). This 
system enabled us to guide and steer amphitrichous Magnetospirillum magneticum strain AMB-1 to 
specific location via magnetic islands. The geometries of the islands allowed us to have control over the 
specific magnetic field gradients on the bacteria. We estimate that magnetotactic bacteria located less than 
1 fim from the edge of a diamond shaped island experience a maximum force of approximately 34 pN, 
which engages the bacteria without trapping them. Our system could be useful for a variety of applications 
including magnetic fabrication, self-assembly, and probing the sensing apparatus of magnetotactic bacteria. 



Controlling magnetic field gradients while simultaneously imaging biological and non -biological magnetic 
specimens will contribute to scientific fields ranging from magnetic fabrication^'^ to probing the behavior 
of magnetotactic bacteria^. A number of techniques currently exist to create and manipulate magnetic 
fields including magnetic tweezers^, magnetic needles^, and permanent magnets'^. Each of these systems provides 
distinct ways to apply magnetic fields and exert forces for magnetic manipulation. Magnetic manipulation tools 
join the repertoire of tools, such as optical tweezers and atomic force microscopy (AFM), for making measure- 
ments like the contour and persistent length of DNA molecules and the strength of cellular adhesion forces^'^"^. 

Magnetic techniques are often preferred over optical techniques for force application because magnetic fields 
provide a way to non-invasively manipulate cells without excessive heating. Many researchers have attached 
commercially available functionalized magnetic beads to cells to examine such cellular properties and functions as 
elasticity and mechanotransduction^^"^^. In biology, there has been a growing interest in organisms like magne- 
totactic bacteria (MTB), because they naturally synthesize chains of magnetic particles and this unique feature 
makes them guidable to specific locations^^. A significant drawback of optical tweezers is that the local heating can 
irreversibly denature proteins and damage living cells. Although magnetic fields applied at low frequencies 
(<100 Hz) typically do not generate significant heat in biological samples, heat may be generated in magnetic 
coils themselves. When the magnetic coils are located at a distance from the sample, such as with Helmholtz coils, 
the generated heat can be removed by convective cooling with liquid^^ or gas^^ before it affects the sample. 

In this paper, we present an approach that combines macroscale Helmholtz coils and lithographically fabri- 
cated nanoscale islands composed of permalloy (NisoFeao) to control magnetotactic bacteria as well as magnetic 
particles (Fig. la-c). These nanoscale islands concentrate magnetic fields generated by the Helmholtz coils, and 
the resulting magnetic field gradients apply forces on the magnetotactic bacteria. The system controls applied 
forces and moments and enables the precise positioning and orientation of magnetic particles as well as being able 
to steer motile magnetotactic bacteria. In comparison to techniques like magnetic tensile cytometry (MTC)^^, 
which also uses concentrated magnetic field gradients, the use of nanoscale magnetic islands bypasses the need for 
the bulk pole pieces (i.e., large needles) typically used to generate magnetic field gradients^^ and can create arrays 
of islands for high throughput control of magnetotactic motility. The advantage of our approach is applying 
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Figure 1 | Integrated nanofabricated metallic island and optical microscope for control magnetic fields on magnetic particles and magnetotactic 
bacteria, (a) Schematic of experimental set-up: Helmholtz coils integrated with nanofabricated permalloy (Ni8oFe2o) islands, (b) Custom-built pair of 
Helmholtz coils to fit an inverted microscope for magnetizing the permalloy island. Inset (i) is a front view of the Helmholtz coils and the specimen 
holder. Inset ( ii) is a side view of the Helmholtz coils (c) Characterization of the Helmholtz coils showing magnetic field strength versus applied current to 
the Helmholtz coils. 



physics to biology by building a system that controls magnetotactic 
bacteria. This capability is of great interest to many areas including 
(in the short term) investigating the behaviors and mechanisms of 
magnetotactic bacteria (MTB) and (in the much the longer-term) the 
idea of creating sentinel robotic systems that could be directed with 
drug treatments, etc., in the body through controlling magnetic 
fields. While the latter is in the very long term, it is interesting since 
these organisms are non -pathogenic and naturally respond to mag- 
netic fields. This system described here may impact multiple areas of 
research ranging from directed magnetic self-assembly to the invest- 
igation of magnetically responsive organisms and magnetically 
based drug release systems^^. 

Results 

In order to create the magnetic field concentrators for integration 
with the Helmholtz coils (Fig. 1), we fabricated thin permalloy metal- 
lic islands. Figure 2a is a schematic of the fabrication process for the 
permalloy islands (Fig. 2b; see Material and Methods section for a 
detailed description of fabrication of these islands). To visualize the 
interaction of AMB-1 with these islands using a 63X, 1.4 N.A., oil 
immersion objective in an inverted microscope, the islands and the 
microscope slide thickness could not exceed 170 |im, and when the 
thickness of the metal layer was less than 30 nm, the metal layer was 
optically translucent. This fabrication approach allowed us to visu- 
alize the bacteria and magnetic particles on top of the islands. In 



order to exert higher forces on the beads, we also nanofabricated 
islands four times thicker than the 30 nm islands (120 nm). 

Forces on the magnetic beads and the magnetotactic bacteria. We 

examined our system through simulations (Fig. 3a, b, top panels) and 
experimentally compared these using ferromagnetic particles made 
from uniform polystyrene spheres coated with chromium dioxide 
(4.5 jim in diameter; Spherotech Inc.; see supplementary material for 
the code used to run these simulations). These particles permitted us 
to map the magnetic field gradients and quantify the forces imparted 
on the magnetic particles by these islands. The ferromagnetic 
particles moved to the areas of high magnetic field concentration 
(Fig. 3b, bottom panels. Movies SI and S2). The beads moved to 
the regions predicted by finite element method magnetics mo- 
deling (Fig. 3b, top panels). We obtained magnetic hysteresis 
curves to provide important information such as magnetization 
(from which the magnetic susceptibility was determined) and the 
coercivity for our approach. The magnetic hysteresis of the 
magnetic beads and the permalloy film were acquired using an 
alternative gradient force magnetometer (AGFM) (MicroMag 
Model 2900; Princeton Measurements Inc.) and exhibited coercive 
forces (He) of 567 Oe (in agreement with other studies^^) and 14 Oe, 
respectively (Fig. 3c, d). Our results were in line with the known low 
coercive force of this permalloy^^. This magnetic hysteresis also 
indicated that there was little magnetic remanence remaining in 
the islands once a field was removed. 
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Figure 2 | Fabrication of permalloy islands to make magnetic concentrators, (a) The image reversal fabrication steps implemented to make permalloy 
islands of controlled geometries, (b) Bright field microscope images of different shapes of the permalloy islands fabricated using the image reversal 
method. Scale bar, 10 [im. 



To quantify the forces exerted by the islands on the beads, we 
applied Newton's second law and the balance of magnetic and drag 
forces (see supplementary material for equations). Using in-house 
code written in MATLAB, we determined the magnetic forces 
exerted on the beads by tracking them. Figure 4 shows the magnetic 
force with respect to the distance to the edge or vertex in the vicinity 
of the square and diamond islands. We also varied the thickness of 
the island, and found that increasing the thickness from 30 nm to 
120 nm subsequently increased the magnetic force by approximately 
50% (Fig. 4). Generally, islands with larger areas and volumes exert 
more forces on the beads, which is consistent with theory predicting 
that magnetic force is directly proportional to the volume of the 
magnetic material. Furthermore, we found that the diamond-shaped 
island (which are the same area as the square) with its shape aniso- 
tropy created higher gradients. 

In addition, we calculated the forces on the magnetic bacteria by 
extrapolating data from the magnetic beads forces that we have prev- 
iously calculated. We used a similar approach to calculate the force 
exerted on these magnetotactic bacteria by the magnetized islands. 
To do this, we used the magnetic hysteresis (of the beads and mag- 
netotactic bacteria. Fig. SI) and the forces vs. distance curves. From 
the magnetic hysteresis we calculated the magnetic susceptibility. 
From the force versus distance curves, we determined the forces on 
the beads very close to the islands' edge (see the supplementary 
material for details of the calculations and equations used). 
Assuming that the biogenic magnetic particles (BMP) are 1.9 X 
10"^^ m^ (50 nm in diameter) in cross-sectional area and that the 
chains have a maximum chain length of about 2 jim (a cells is typ- 
ically around 3 |im) we estimate that very close to the edge (approxi- 
mately 3 |im), the magnetic forces on the AMB-1 starts to increase 
exponentially and reach a maximum of 34 pN (Fig. 5). Figure 5 also 



shows the effect of increasing the length of the magnetosome chain. 
For instance, a cell with a chain that is only 0.5 jim would experience 
10 fN (femto Newton) at 1 |im away, however a cell with a magneto- 
some chain that is 2 microns in length would experience 10 fN at 
about 2 jim away. The maximum force on the magnetosome chains 
by the diamond island of sizes 0.5 jam, 1 jim, 1.5 jim and 2 jim is 
approximately 8 pN, 16 pN, 25 pN and 34 pN, respectively. 

Guiding the cells using the notch and diamond shape islands. We 

used our magnetic field concentrators to affect the motion of living 
organisms as well. Since magnetotactic bacteria sense and respond to 
the magnetic field of the Earth, we applied a localized magnetic field 
gradient to motile AMB-1. While magnetic particles passively move 
to locations of the steepest magnetic field gradients, the bacteria 
actively swim in response to the magnetic field gradients. One 
interesting observation was the deflecting behavior (Fig. 6, left 
panel), which we attributed to the "bridging" of the magnetic field 
lines across the "notch- shaped" islands as shown in Figures 6a-b, 
Figure S2 and Movie S3. In Figures 6c, d, and e, the deflection of the 
AMB-1 away from and near the notched part of the island revealed 
distinct deflection responses when using a very low external field 
(3 Oe). Near the notch region, this deflection varies from 0 to 
approximately 50° from the x-axis (Fig. 6d, black symbols) 
whereas in the upstream region the AMB-1 cells tend to orient to a 
lesser degree (0-30°). (Fig. 6e, black symbols). These data suggest 
that near the notch region the cells are deflecting their motility paths 
to follow the magnetic field lines. These deflection patterns are 
different than that of the control (when we apply no external 
magnetic field) as the cells are showing wider spread of deflections 
in both upstream and near the notch regions (Fig. 6d and e, blue 
symbols). These results are important because the method could one 
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Figure 3 | Characterization of the permalloy island and magnetic bead response, (a) A two-dimensional field intensity profile of the application of 
magnetic fields on the permalloy islands through our Helmholtz coils simulated using finite element method magnetics (FEMM). The red circle 
represents the approximate area of homogenous magnetic field where the permalloy islands are placed in the simulation, (b) Simulation results of 
the islands with appUed external magnetic fields from the Helmholtz coils of about 80 Oe (top). Experimental comparison of these concentrated magnetic 
field regions through tracking ferromagnetic bead movement with the applied magnetic field (bottom). Scale bar, 10 |am. Magnetic hysteresis loop 
of the ferromagnetic beads and permalloy film using an alternative gradient force magnetometer (AGFM). (c) AGFM hysteresis loop of the ferromagetic 
beads with a coercive field of 567 Oe and (d) AGM hysteresis loop of the permalloy films with a coercive field of 14 Oe. Both magnetic hysteresis 
curves are normalized by the moment. 



day supplant conventional magnetic particle inspection techniques 
(MPI)'l 

Moreover, we were able to control the movement of the AMB-1 
using the diamond island. The cells are able to move along the edge of 
this island, following it, until reaching a point of a low field or an area 
of a corner where the bacteria then moves away from the diamond 
(Fig. 6, right panel). The island used here is 120 nm thick, and the 
external field strength is 120 Oe. For example, in Figure 6 h-m 
(Movie 84) when the external magnetic field is on, a MTB gets direc- 
ted toward one of the magnetic island vertices labeled as "Engaging 
zone" (this is where the field is the highest based on our simulations; 
Fig. 6f). This MTB is motile, and so it is actively moving, but it is 
moving in a directed fashion. The bacterium appears to follow the 
edge of the island until it reaches another vertex labeled as "Escape 
zone" (the field is lowest in this area based on our simulations; 
Fig. 6f). Thus, we were able to guide the cells (at an angle from the 
external magnetic field) to specific points near the islands, which is a 
significant capability for directing the bacteria. 

To further quantify these effects, we analyzed the videos and deter- 
mined that 30% ± 10 (s. d.) (n = 83) of the cells swimming near these 
diamonds island (we examined a 60 |im by 40 |im region enclosing 
the diamond island) effectively altered their direction toward the 



edge of the islands. The cells that were pulled towards to the edge 
were the ones that were the closest to the edge. We also quantified the 
distance between the cells and the islands when the cells were influ- 
enced by the magnetic field (the moment at which they lose inde- 
pendence in swimming), which was 3.6 jim ±1.4 (s. d.) 

Discussion 

We used our microscope nanofabrication approach to control the 
motility of magnetotactic bacteria (Fig. 6). This control was enabled 
by implementing a combination of a single pair of Helmholtz coils 
and nanofabricating magnetic permalloy islands. Here the external 
magnetic field produced by the coils influences the cells to move from 
left to right or right to left within the field of view of the microscope. 
Thus the external field magnetically induces their orientation. The 
direction of movement of MTB is also controlled when using static 
Helmholtz coils and different geometries of the magnetic permalloy 
islands. These provide local control through local field gradients. 
These movements are influenced by the re-orientation due to the 
localized magnetic fields coupled with the velocity of the flagella 
enabled swimming bacteria. The distance that the cells needed to 
engage with the edge of the diamond island was below 3.6 jam. 
This distance coincides with the rapid increase of the field (Fig. 5). 
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Figure 4 | Force versus distance results for the square- and diamond- shaped magnetic islands with an applied magnetic field. The square and 
circles represent the data for fields equal to 40 G and 120 G, respectively, (a-b) Force due to square island with height equal to 30 nm and 120 nm. 
(c-d) Force due to the diamond island with height equal to 30 nm and 120 nm. 



To increase the area or volume near the edge that effectively engages 
the cells, the system could be influenced by changing parameters 
such as: increasing the current to the coils from the power supply, 
increasing the volume (either the area of the thickness) of the island, 
and by having the cells being more magnetically responsive. This 
could be accomplished through connecting power supplies in series 
or growing the cells under strict anaerobic condition (e.g. by degauss- 
ing their culturing medium). 

This system could also be used to calculate the forces exerted by 
the flagella by considering three limiting cases. If the M. magneticum 
strain AMB-1 are swimming away from the island, and at a critical 
distance from the island, AMB-1 can no longer resist the magnetic 
force then Ffiageiia ^ Fmag- Closer to the island than the critical dis- 
tance, Fjnag will be greater than Ffiageiiaj and the magnetic force will 
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Figure 5 | Forces on the magnetosomes from the islands and the effect on 
increasing the magnetosome chains size, (a) The force vs. distance curve 
distance curves indicate the different sizes in the magnetosome chain of 0.5 
(blue), 1 (cyan), 1.5 (magenta) and 2 (green) microns in length. The 
surface area is fixed and is 1.9 X 10~^^ (50 nm in diameter). The forces on 
the magnetosome increase exponentially within a few microns from the 
edge of the island. 



pull microorganism towards the island. The third limiting case 
occurs when the bacteria overcomes the magnetic forces exerted by 
the islands (F^ag < Fflageiia)- However, this calculation could only be 
made in an ideal case scenario, and it is more complicated than 
looking at the limiting cases because the cells may be slowing down 
due to environmental cues as they response to many difference sig- 
nals. The force on the AMB-1 about 2 |im from the edge of the 
diamond is about 10 fN (assuming the length of the chain is 2 jim 
long). The estimated maximum force of 34 pN is in indirect agree- 
ment with the estimate force on the flagella motor of other species 
{Escherichia coli and Salmonella Typhimurium) which is about 
1500 pN^^. This force of 34 pN is sufficient to engage the cells, but 
not to complete immobilized them thus we can effectively guide these 
cells. 

One important potential future concept from these additional 
studies is that these islands could be used in single cell manipulation 
to screen MTB with weak and strong magnetic chains by sorting 
them with an array of diamond shape islands position diagonally 
from each other such as to a compartment in a microfluidic channel 
device. Such arrays of islands could also be potentially used in direc- 
ted evolution by trapping the cells repeatedly. 

To conclude, we developed a custom-designed system integrating 
Helmholtz coils and nanofabricated magnetic islands for applying 
controlled magnetic fields and field gradients. We used this system to 
control the movement of magnetic particles and to study the beha- 
vior of magnetotactic bacteria. Straightforward modifications of this 
system are possible to exert higher forces and would be useful for 
studying cells and other organisms that sense magnetic fields. Our 
work has implications for several fields including magnetic self- 
assembly, physiological response to magnetic fields, and research 
on organisms containing magnetic particles. 

We would like to thank Carson Kline for his initial help in the 
NanoFab. The authors would also like to thank Andrew Gamble for 
help with measuring the hysteresis loops and Dr. Sara Majetich for 
helpful comments on the manuscript. This work was supported in 
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Figure 6 | MTB interactions with notch and diamond islands. (Left panel) MTB were tracked as they moved from the upstream region to the notch 
region, (a) Figure showing the different regions (upstream and notched regions), (b) Schematic showing the magnetic field lines near the notch island, 
(c) A schematic of the coordinates to track the deflection. Orientation of the MTB when they are upstream of the notch, Uj. (d) and as they approach the 
region above the notch, Uj- (e) when compared to vertical distance (y) away from the islands. The field corresponding to d-e is 3 Oe and the thickness of 
the island is 30 nm. (Right panel). The black symbols represent the bacteria when a field is applied and the blue symbols represent the bacteria when there 
is no external magnetic field. (Right panel) This panel shows the movement of MTB along the edge of the diamond island, (f) FEMM simulation 
of the diamond magnetic island showing the differences in the field between the "engaging" and "escape" zones, (g) A micrograph of the diamond island, 
which is 120 |am in thickness and the external magnetic field is 120 Oe. (Scale bar, 5 |am). (h-m) Image sequences of the magnetotactic bacteria moving 
along the edge of the diamond island. The bacterium (highUghted with a red circle) moves toward the engaging zone and gUdes along the edge 
(as in a "rail") of this island until moving away from the island at the escape zone. 
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Methods 

Construction and characterization of Helmholtz coils. The Helmholtz coils in our 
system produce a uniform magnetic field with field lines concentrated within the 
permalloy islands. The governing equation used to calculate the forces on the 
magnetic particles is equation 1 although it can be found in other forms with the 
assumption that the initial magnetization of the magnetic beads is negligible^^'^*^. 

F=(m-V)B (1) 

To develop an approach for manipulating magnetic samples including magnetic 
particles and magnetotactic bacteria, we first custom designed a pair of Helmholtz 
coils for a conventional inverted microscope (Zeiss Anxiovert 200) as shown in 
Figure 1. The axis of these coils was positioned perpendicular to the optical axis. This 
pair of Helmholtz coils consisted of two identical coils separated by the radius of the 
coils. We designed and fabricated a custom stage that held the samples and the 
Helmholtz coils through an L- shaped piece attached to a micromanipulator 
(Newport, M-462) (Fig. lb). The frames of the Helmholtz coils and stage were made 
of poly ( methyl- methacrylate) (PMMA). We used a pair of miniature cooling fans to 
prevent excessive Joule-heating. The radius of the coils were 4.5", and the gauge of the 
magnet wire used was 22 AWG (American wire gauge). This pair of Helmholtz coils 
had 400 turns in each coil, and the windings were stacked in 7 layers. The magnetic 
field was measured using a gaussmeter (Lakeshore, Model 410). To characterize the 
magnetic field, B, at the center location of the coils, we used 
equation 2: 



where N is the number of turns, |j,o is the permeability constant, I is the applied 
current, and R is the radius of the coils. Figure Ic depicts the characterization of the 
field with respect to the current. 

Fabrication of the custom-made metallic magnetic islands for an inverted 
microscope. A 150 |im glass wafer (Mark Optics Inc.) served as the microscope slide. 
The image reversal technique resulted in a negative image of the mask pattern, 
allowing for effective lift-off of the metal film after fabrication^^. Using heat-releasing 
double-sided tape (release point = 150°C), a thin, fragile glass wafer was attached to a 
stronger silicon wafer for generating the islands. Bis(trimethylsilyl)amine (HMDS) 
was spun onto the glass/silicon wafer sandwich twice at 4000 rpm for 45 s. HMDS is 
an adhesive promoter commonly used to adhere photoresist to a glass wafer. A 
positive photoresist, AZ5214E-IR (MicroChem), was spun onto the surface at 
4000 rpm for 45 s. The wafer was soft-baked at 95°C for 60 s to dry the solvent. 
During the initial exposure, the resist was controllably exposed to ultra-violet (UV) 
light through a photomask with an intensity of 150 mj/cm^ for 50 s. The photomask 
was a chrome quartz mask with a tolerance of 0.5 \im, which enabled us to repeatedly 
fabricate sharp features. To activate the crosslinking agent exposed to the UV light, 
the wafer was incubated at 105°C and 115°C on hotplates for 30 s and 120 s, 
respectively. During the flood exposure, the wafer was exposed to UV light without 
the mask at 150 mj/cm^ for 150 s. This step makes the photoresist not exposed in the 
first round of UV light soluble when introduced into the developer. The resist was 
then developed with a developing solution ( AZ400K, Clariant Corporation) diluted in 
deionized (DI) water (1:4) for 0.8 min followed by lightly rinsing with DI water and 
drying with nitrogen. 

The ensuing master was used as the mold to deposit the metals (Fig. 2a, metal 
deposition step). Using a sputter deposition system (Consolidated Vacuum 
Corporation), a 25 A layer of tantalum and a 250 A thick layer of permalloy were 
subsequently deposited onto the master mold. The tantalum layer helped with 
bonding of the permalloy to the glass wafer. The glass/silicon sandwich was then cut 
along cutting lines defined in the photomask to produce coverglass-like pieces that fit 
into a custom-made coverglass holder. The cut, coverglass- sized silicon pieces with 
deposited islands were released from the silicon wafer upon heating the wafer to 
150°C. Excess photoresist and permalloy were removed by immersing the samples in 
acetone until the film residue completely lifted off from the glass substrate 
(<15 min). As a final step, the glass wafers were rinsed in fresh acetone, isopropyl 
alcohol, and DI water. Figure 2b shows brightfield images of different fabricated 
islands used in this study. 

Microscopy and imaging. We mounted the Helmholtz coils on an inverted optical 
microscope (Zeiss Axiovert 200) equipped with a halogen lamp and an oil immersion 
objective (Zeiss 63 X, N.A. = 1.4). To initially find the island we used a 5X objective 
and then a 63X oil immersion objective. To track single magnetic beads, we diluted 
the magnetic particles from the manufacturing suspension (4.5 |j,m in diameter; 
Spherotech Inc.) with DI water to get 1% v/v. We captured the images using a CCD 
camera (Insight QE, SPOT Imaging Solutions) and analyzed the images using ImageJ 
and MATLAB Image processing tools. 



Magnetotactic bacteria cultivation. Clones of Magnetospirillum magneticum strain 
AMB-1 (Magnetospirillum sp. ATCC® 700264) were obtained by aseptically 
streaking an incoculant across an agar plate containing magnetic spirillum growth 
media (MSGM) with the same nutrient components as the liquid media. Individual 
dark colonies were selected and microaerobically grown on liquid media until the 
beginning of the stationary phase. The liquid media contained 10.0 mL of Wolfe's 
Vitamin solution, 5.0 mL of Wolfe's Mineral solution, 2.0 mL of 0.01 ferric quinate, 
0.45 mL of 0.1% resazurin, and various salts. These salts included 0.68 g of KH2PO4, 
0.12 g of NaNOs, 0.035 g of ascorbic acid, 0.37 g of tartaric acid, 0.37 g of succinic 
acid, and 0.05 g of sodium acetate. To prepare the 0.01 M ferric quinate, we used 
0.27 g of FeCls, 0.19 g of quinic acid, and 100.0 mL of deionized water (DI) water. 
Resazurin was used as an oxidation-reduction indicator to visually monitor the level 
of oxygen present in the liquid media. The fading of the initial pink color of the liquid 
media indicated that the bacteria were in exponential growth phase. To confirm 
magnetic particle uptake within growing bacteria, we visually inpected the cells 
following placement of the screwcap culture tubes on a magnetic stirrer and observing 
the scattering of light, as the ambient light is reflected off these particles. We also 
inspected the AMB-1 with a transmission electron microscope (TEM). 

Finite Element Methods Magnetics (FEMM). We ran the simulation using freeware 
- Finite Element Methods Magnetics (FEMM). We set up an asymmetrical 2D 
problem by simulating half of the Helmholtz coils as shown in 3A. We wrote an in- 
house Lua script (a programing language that can be interpreted by the FEMM 
program) to run simulations with different geometries (square, diamond and notch). 
We fixed the radius of the coils, the number of turns in the wire and the position of the 
islands and varied the current in the wire (see details of the script in the 
supplementary material). 
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